Generalized Chiral Symmetry and Stability of Zero Modes 
for Tilted Dirac Cones 
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While it has been well-known that the chirality is an important symmetry for Dirac-fermion 
systems that gives rise to the zero-mode Landau level in graphene, here we explore whether this 
notion can be extended to tilted Dirac cones as encountered in organic metals. We have found 
that there exists a "generalized chiral symmetry" that encompasses the tilted Dirac cones, where 
a generalized chiral operator 7, satisfying 7^_ff -I- = for the Hamiltonian H, protects the 
zero mode. We can use this to show that the n = Q Landau level is delta-function-like (with no 
broadening) by extending the Aharonov-Casher argument. We have numerically confirmed that a 
lattice model that possesses the generalized chirality has an anomalously sharp Landau level for 
spatially correlated randomness. 
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Introduction — Since the seminal observation of the 
quantum Hall effect in graphene, [H 0] the zero modes 
of the massless Dirac fermion, which is essential for the 
characteristic quantum Hall effect in graphene, has been 
intensively discussed. Specifically, the zero mode is sta- 
ble against ripples in graphene, which has been discussed 
in terms of the index theorem, or more explicitly for 
wave functions with the argument due to Aharonov and 
Chasher For the stability of zero modes, a cru- 

cial ingredient is the chiral symmetry defined as the ex- 
istence of an operator F that anti-commutes with the 
effective Hamiltonian H , {F, H} = 0, with F^ = 1, which 
is more relevant than the Dirac-cone dispersion per se, 
since the stability is also shown for a chiral symmetric 
bilayer graphene with a quadratic dispersion [lo| . In 
the case of graphene, the Dirac cone is vertical, so that 
the chiral operator has a simplest possible form of F = cr^ . 
Even when the system is disordered, the anomalous crit- 
icality at the n — Landau level is retained as far as the 
randomness respects the chiral symmetry (TT| - [l^ . The 
stability of zero modes has been observed experimentally 
as well for a mono-layer graphene 
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On the other hand, we encounter Dirac cones in a wider 
class of materials, where the effective theory is more gen- 
erally described by tilted Dirac cones, as is the case with 
an organic material a-(BEDT-TTF)2l3 [H-il. Since 
the conventional chiral symmetry is broken in tilted Dirac 
cones, it becomes an essential question to ask whether 
(i) the symmetry is entirely broken, and (ii) whether the 
anomaly in the systems with usual chiral symmetry is 
washed out in tilted Dirac cones. In the absence of dis- 
order, the eigenvalues and eigenfunctions for tilted Dirac 
cones in magnetic fields have already been obtained in 
existing literatures (23l - |25j , which indicate that the zero- 
energy modes themselves persist in clean tilted cones. 

Effects of disorder on the zero modes is indeed an im- 
portant issue, since, for graphene with the conventional 
chiral symmtry, it has been established that the disor- 
der has an anomalous effect of retaining a sharp (delta- 



function-like) zero Landau level accompanied by a sharp 
(step- function-like) quantum Hall step, if the disorder re- 
spects the chiral symmetry. We can thus pose a question: 
what is the effect of disorder for the zero Landau level 
in tilted Dirac cones, where the conventional chiral sym- 
metry is absent. This is exactly the motivation of the 
present work. In particular, we shall look at how the sta- 
bility of the zero modes possessing an anomalous critical- 
ity at the n — Landau level in the presence of disorder 
(e.g., random components in magnetic fields) is affected 
by the breakdown of the conventional chiral symmetry in 
tilted Dirac cones. Curiously, the study leads us to find 
a "generalized chiral symmetry" , which is then shown 
to give rise to an even wider stability of the zero modes 
persisting in tilted Dirac cones. 

Thus we shall first generalize the conventional chiral 
symmetry so that the tilted and untilted Dirac disper- 
sions can be captured in a unified manner. The gener- 
alized chiral symmetry is defined by the existence of an 
operator 7, which is not necessarily Hermitian, satisfying 
the relation 'y^Hj = —H and 7^ = 1. This is consistent 
with the condition for the Dirac operator to be elliptic 
that is required for the index theorem [26*]. It enables 
us to show a topological stability of zero modes that is 
present in tilted Dirac cones. With this generalized chi- 
ral symmetry, we reformulate the eigenvalue problem so 
that the Aharonov-Casher argument for counting the 
number of zero modes is extended to tilted cones. This 
implies the zero Landau level is indeed delta-function- 
like. We then numerically confirm how these field theo- 
retic treatments on the stability of zero modes appears in 
a lattice fermion model that has tilted Dirac cones where 
the titling is varied continuously. 

Formalism — The effective Hamiltonian for a tilted 
Dirac cone can be generically expressed as [23l - [25| 



H = (7QiW ■ 7T/h) + (a- ■ X)7r^/h + (cr • Y)Try/h, (1) 

where (cTa:, f y, c^) are Pauli matrices while ctq is a 2 x 2 
unit matrix, and X,Y , W are real. The Dirac cone is 
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tilted when *W = {Wx,Wy) is nonzero, while *X = 
{Xx,Xy,Xz) and *Y — {Y^,Yy,Yz) characterize the 
anisotropy of the Dirac cone. Here tv = p + eA (e > 0) 
is the dynamical momentum with *p = —ih{dx,dy) 
and the vector potential *A — {Ax,Ay) for the mag- 
netic field B = dxAy — dyAx perpendicular to the x- 
y plane. Note that the dynamical momentum satis- 
fies the commutation relation [TTa^jTTj,] = —iheB. For 
the conventional Dirac cone {W = 0), the Hamilto- 
nian has the chiral symmetry associated with an oper- 
ator defined (and generalized to anisotropic cases) by 
r = cr • {X X Y)/\X X Y\, which anti-commutes with 
the Hamiltonian H = [cr ■ X)Trx/h + (cr ■ Y)'Ky/h with 

When the Dirac cone is tilted (with H containing ctq 
for W 7^ 0), the conventional chiral symmetry is broken. 
However, here we find that a generalized symmetry does 
exist, which can be seen if we introduce a new operator 
7 defined by 

7 = cr • [(X X r) - i{WyX - WxY)]/A 



with A2 = |X X Y\^ - [WyX - WxY)"^. This operator 7 
is non-Hermitian for W ^ 0, but its eigenvalues are ±1 
since 7^ = (7^)^ = (Jq. We can show that, if > 
satisfies a relation with the Hamiltonian, 
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which we call the generalized chiral symmetry. This 
symmetry reduces to the conventional chiral symmetry 
(7 — >■ r) for W = 0. It should be noted that the cross 
section of the tilted Dirac cone with a constant energy 
plane is an ellipse as long as > (while a hyperbola 
when A^ < 0). 

Can we say anything about the wave functions as a 
direct consequence of this generalized chiral symmetry? 
For this purpose it is instructive to choose the (right-) 
eigenvectors |±) of the opeator 7 (with 7|±) = ±|±)) 
in the spinor space as a basis of the eigenvalue problem 
for H . If we express the (normalized) wave function as 
-0 = + |— )?/'_, the Schrodinger equation Hip = Etp 

reduces to 
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Note that the operator 7, being non-Hermitian, has 
eigenvectors that are in general not orthogonal with 
each other with j3 = (+|— ) 7^ 0. When the gener- 
alized chiral symmetry, "f^ = —H, holds, we have 
{+\H\+) = {-\H\-) = 0. The Schrodinger equation 
then becomes 
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where we have introduced a complex *« 
h-'{{+\Wxao + X ■ {+\Wyao + Y 



= [ax, ay) = 
cr|-)). The 



Schrodinger equation for the E = states (zero modes) 
therefore amounts to that for the zero modes of the un- 
tilted (but can be anisotropic) Dirac cones (with zero 
diagonal elements). We see that the zero mode is then 
given by 



a ■ Tvip- — and a* ■ Tvip+ = 0. 



(2) 



As we shall see in the following, the chirality ± has to 
be assigned to (normalizable) wave functions iIj±. We 
can also note that the equations for zero modes apply to 
spatially varying magnetic fields B{x,y) as well, which 
can even be random. 

Aharonov-Casher argument extended to the general 
chirality — Is the stability of zero modes inherited by 
the general chiral symmetry? For this we can look at 
the zero modes for the tilted Dirac cone in a spatially 
varying magnetic field B(x,y), which are therefore de- 
termined by eq.©. Here it is convenient to adopt "prin- 
cipal coordinates" by rotating a with an orthogonal ma- 
trix T (with det r = 1), so that the complex numbers 
(zxjZy) — {ax,c!tY)T~^ become orthogonal with each 
other on the complex plane. It is indeed possible to do 
this with {zy/zx)/\zy/zx\ = -sgn(Im(Q!jf ay))z = ~ix- 
With the transformed dynamical momentum 11 = Ttv 
and the coordinates R =* {X,Y) ^ Tx, Eq.© for -0_ 
reads 



a ■ Tv-ijj^ — Zx\ ""^(AH 
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where the " ellipticity" A = \/\zx\/\zy\ > is a positive 
real number. The equation for ip- then becomes (AHx — 
«xA~^nY')'0_ — 0. This, along with a similar equation 
for ■(/;+, reduces to 



(AHx±zxA-iHy)V± = 0. 



(3) 



With these equations, it is straightforward to general- 
ize the analytic Aharonov-Casher argument (sj for the 
zero modes to the present case. We use the "Coulomb 
gauge", X^dxAx + A~^9yAY = 0, which is automati- 
cally satisfied if we introduce a "scalar potential" (p with 
(Ax, Ay) = i~\-^dYip,X^dx(p)- Then Eq. © simpli- 
fies to 



ih 



V^± =0 



with 2?± = (A9x ± ixA ^9y) and (j>o — h/e the flux 
quantum. Putting ^± — exp(±27rx(^/(/)o)'0±, we finally 
obtain 

with R = iX,Y) = {X-'^X,XY). Namely, the function 
■0-1- is an entire function of Z± — X ± ixY over the whole 
complex plane, namely a polynomial in Z±. The func- 
tion ip is determined by B^ — BxAy — OyAx = {^^d\ + 
A^^c)y)V3 = (9^ -I- i9|,)i^, which implies that ip(R) = 
J dRG{R- R')B,{R') with G{R) = (l/27r) log(r/ro) 
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and — B? . When the magnetic field is nonzero 
only in a finite region, we have an asymptotic behavior, 
($/27r) log(r/ro), for r — > oo, where $ = / dRB^ — 
J dR Bz is the total flux. Then we obtain 

for r — >■ cx), in which ±x$ < is necessary for ^± to 
be normalizable. The normalizability of |V'±lf, indicates 
that the degeneracy of zero modes is ^/^^o 0- This is 
exactly equal to the total number of energy levels in the 
Landau level, which implies a remarkable property that 
no broadening occurs for the n = Landau level. 

Higher Landau levels — We can also note that the 
present formulation provides a simple algebraic repre- 
sentation of the Landau levels, including higher ones, 
in a uniform magnetic field B{x,y) = B > 0, where 
the Schrodinger equation reads a ■ ^ip^ ~ Ei-'+ and 
a* ■ 7rip+ = Ei/j-. Here tt = tt — q with a real q sat- 
isfying a ■ q = E/3, which amounts to taking the origin 
to be the center of the cross section (an ellipse) of the 
tilted Dirac cone on a constant-energy plane. Since we 
have [Tr^^jTTy] = [712^,715^], an annihilation operator cle sat- 
isfying [a£;,a^] = 1 can be defined as a^; = a ■ tz/Ab 
with Ab = y/2lm{axoip)heB when x = 1- With a 
basis fi^E — {a-EYfo/y^- for positive integers £ where 
acting as the usual raising operator with ciEfo = 0, 
the eigenstate Hip„ = E^ipn with a nonzero energy 
En = sgn(n)-\/|n|AB, found in Refs. p3l425i |. can be 
algebraically expressed as — sign(n)|-|-)/|„[_i^^^ -|- 
|-)/|n|,ii;„ forn = ±l,±2,.... 

A lattice model — In actual materials with general 
band structures, the above argument based on the effec- 
tive Hamiltonian ([T]) and the generalized chiral operator 
7 holds only as a low-energy effective model. In order 
to confirm whether the anomalously sharp zero-Landau 
level discussed above appears as well in a lattice model 
that possesses tilted Dirac cones at low energies, we con- 
sider a two-dimensional lattice model with a Hamilto- 
nian having nearest (t) and second neighbor (t') hoppings 
(FigHI) on a square lattice, 

^ = E + h.C.) + {-ir+n{cU,^Cr + h.C.) 

r 

+ E *'(4+e.+e„Cr + 4+e.-e,Cr + h.c), 
r 

where r = (x, y) denotes a lattice point in units of the 
lattice constant, and e.j:{ey) the unit vector in x{y)- 
direction. The primitive vectors for the present lattice 
system can be chosen as ei = — Ry and e2 — Ex -'r Cj, 
(FiglT]). In the absence of magnetic fields, the Hamilto- 
nian in the momentum space becomes 

TT/,\_ 2<'(cosA:i -I- cosfc2) ^(fe) 

~ |_ A*(fc) 2t'(cosfci -Hcos/fcs) 

with A(fc) = -t(-l + e^*-^i + e*'=i+*'=2 -I- e*'^^), where 
ki = k ■ ei and k2 — fe • 62. The band dispersion 



in the first Brillouin zone has a pair of Dirac cones at 
{ki,k2) = (7r/2,-7r/2) and (-7r/2,7r/2) with £; = as 
long as \t'/t\ < 0.5, where the tilting becomes stronger 
with \t'\ (Fig. [1]). We can then show that the effective 
Hamiltonian around these Dirac cones is reduced to the 
Hamiltonian ^ with *W = (0,±4t'), *X = (0,2t,0) 
and *Y ~ (=F2t, 0, 0). The model reduces to the 7r-flux 
modcl[27] when t' = (FiglUb)). 




by the primitive vectors ei and 62. Energy dispersions E{k)/t 
for (b) t' = 0, (c) t'/t = 0.2 and (d) t' /t = 0.4. 

We apply a magnetic field to this model to examine the 
stability of zero modes against the disorder in the mag- 
netic field. The magnetic field is taken into account by 
the Peierls substitution t te-'^"^'^'^\ t' f f^-^T^^e' [r] 
such that the summation of the phases 0(r), 9'{r) around 
a loop is equal to the encircled magnetic flux in units 
of the flux quantum. Here we have adopted the string 
gauge [28J to treat smaller magnetic fields. Disorder is 
introduced here as a random component, S(j){r), in the 
magnetic flux 0(r) — (f> + S(t>(r) piercing each plaque- 
tte, where (f) is the uniform part. The random part 
S(j>{r) is assumed to have a gaussian distribution with 
a variance a and a spatial correlation length 77 with 
((50(r)(50(r')) = (<502)exp(-|r--r'|V477) Hg. We have 
chosen this disorder since it restores, for large enough 77, 
the generalized chiral symmetry of the effective Hamilto- 
nian at tilted Dirac cones. 

Figure[2]displays the density of states for t'/t = 0.4, for 
which the tilting is signiflcant (FigUId)). The result is 
obtained by the exact diagonalization of a flnite system, 
with an average over 5000 samples is performed. It is 
clearly seen that the n — Landau level becomes anoma- 
lously sharp even in random magnetic fields as soon as 
the correlation length 77 of the random flux exceeds the 
lattice constant a, while other Landau levels are broad- 
ened in a usual fashion. This anomaly, appearing only for 
the n = Landau level, suggests that the n = Landau 
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FIG. 2. (Color Online) Density of states for the model with 
t' jt = 0.4 depicted in the inset is plotted against the spatial 
correlation length of the random component of the magnetic 
field 77 for a uniform magnetic field 0/00 = 1/100 and the 
amplitude of the random magnetic field a/^o = 0.0029. The 
result is an average over 5 x lO"' samples with the system-size 
30a by 30o. 



states are degenerated at i? = 0, endorses the stability of 
zero modes for tilted Dirac cones. 

In summary, we have found that the conventional chi- 



ral symmetry can be extended to a generalized chiral 
symmetry that encompasses the models having tilted 
Dirac cone, so that the untilted and tilted Dirac cones 
can be treated in a unified way. The stability of zero 
modes can be proved under this generalized chiral sym- 
metry with an Aharonov-Casher argument. We have fur- 
ther shown, numerically for a lattice model, that topo- 
logically protected zero modes of tilted Dirac ferniions 
survive even in random magnetic fields correlated over 
a few lattice constants. These results suggest that the 
anomaly at n = Landau level can be observed gener- 
ally in systems with tilted Dirac dispersions. 

As a significance of this, we can finally note that the 
existence of the generalized chiral symmetry (A^ > 0) 
is equivalent to the ellipticity of the Hamiltonian ([T]) as 
a differential operator, under which the index theorem 
[26j can be applied. It is an interesting future problem 
to extend the notion of the generalized chiral symmetry 
to a broader class of Dirac-cone systems [soj . 
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